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We consider single-particle quantum transport on parametrized complex networks. Based on 
general arguments regarding the spectrum of the corresponding Hamiltonian, we derive bounds for 
a measure of the global transport efficiency defined by the time-averaged return probability. For 
tree-like networks, we show analytically that a transition from efficient to inefficient transport occurs 
depending on the (average) functionality of the nodes of the network. In the infinite system size 
limit, this transition can be characterized by an exponent which is universal for all tree-like networks. 
Our findings are corroborated by analytic results for specific deterministic networks, dendrimers and 
Viscek fractals, and by Monte Carlo simulations of iteratively built scale-free trees. 
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Introduction - Complex networks are intensively used 
as models for a panoply of physical, chemical, biological, 
or sociological systems, see [1-3] and references therein. 
They have been proven to be extremely useful in under¬ 
standing the statistical as well as the dynamcial features 
of these systems, e.g., the complexity of the structure of 
the internet is captured by scale-free networks, where the 
probability distribution of the number of connections of 
a given node (hub) follows a power law [4]. In contrast 
to stochastic networks, also deterministic networks - such 
as regular, hyper-branched, or fractal types - have been 
applied to study the properties of, say, macromolecular 
compounds and polymers [5, 6]. 

Interestingly, many topologically different networks 
show similar statistical (equilibrium) properties, i.e., they 
fall into the same universality class [3]. Recently, it 
has been shown that, in addition to these classes, there 
are also dynamic universality classes, referring to simi¬ 
lar classical diffusive dynamics’ of different type and on 
different networks [7]. However, it is by no means clear 
that a corresponding quantum dynamics on such net¬ 
works allows for a classification based on a quantum ana¬ 
log of dynamic universality. As will be presented below, 
one can define such a quantum analog also for (single¬ 
particle) quantum transport processes on complex net¬ 
works. In particular, in this paper we will be concerned 
with tree-like networks, i.e., with networks without loops. 
For these we find that one can indeed define classes 
which solely depend on the spectrum of the correspond¬ 
ing Hamiltonian. As we show, classes of parametrized 
Hamiltonians lead to a power-law dependence of a (time- 
averaged) characteristic quantity describing the trans¬ 
port efficiency [8]. 

Quantum Transport Efficiency - We consider closed, 
microscopic physical system modelled by undirected 
graphs/networks Q(V,£) having a set V of vertices/nodes 
and a set £ of edges/bonds. To each node j, out of 
the total number of N, we associate a (basis) state |j), 
such that all these states form a basis set of the Hilbert 
space. The direct couplings between two nodes j and 


k are mediated by a bond; the number of direct cou¬ 
plings of node j is called functionality (or degree) fj. 
The dynamics starting from a localized state \j) is then 
modelled by Schrodinger’s equation [h = I) for the tran- 
siton amplitudes ak,j{t) = {k\exp(—iHt)\j), which is 
properly solved by diagonalizing H [8, 9]. We will con¬ 
sider only such Hamiltonian with indentical coupling 
strengths j = {k\ H |j) = I between any pair of nodes 
connected by a single bond and with on-site potentials 
Hjj = H{fj), i.e., nodes of the same functionality have 
the same potential. Matrices typically associated with 
graphs Q, like its connectivity matrix = fj) and 

its adjacency matrix {H{fj) = 0 for all j), are contained 
in this class. 

Based on the transition amplitudes ak,j(t) and the cor¬ 
responding transition probabilities TTkjit) = |afcj(t)p, 
we will characterize a network’s transport efficiency by 
the average return values d{t) = ^ 

^(t) = P]- Generally, values of Tt{t) = 

0{1) for almost all t imply - on average - a high prob¬ 
ability for an excitation to remain at the initial node, 
thus indicate inefficient transport, values of 7r(t) ^ I 
for almost all t suggesting the opposite. Clearly, one 
needs the entire knowledge of H’s eigenspace. However, 
for 7f(t) this can be circumvented by using the Cauchy- 
SCHWARZ inequality to obtain a lower bound [8] 

I d{t) p(A)p(£;')e-(^-^')‘ < 7r{t) , (1) 

E,E' 

where \d{t) solely depends on the (discrete) spectral 
density g{E) of H. An asymptotic time-independent 
measure for the global transport efficiency can be intro¬ 
duced by the return-quantities’s infinite time limit 

X= lim - [ dt'\d{t')f < lim - [ dt'Tr{t') . (2) 

t->oo t Jq t->oo t Jq 

As a lower bound, x is most instructive if x = G>(1), 
thus the RHS of Eq. (2) is also of order 0{1). Then we 
regard the global transport as being inefficient, whereas 
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for X ^ C>(1) there is no strict implication on the exact 
value of the infinite time limit of Tt(t). However, previous 
results suggest that in particular the maxima of 7r(t) are 
well reproduced by the lower bound |d(t)p, therefore, 
also indicating that the values of the RHS of Eq. (2) lie 
close to the values of x [8] • 

Combining Eqs. (1) and (2) allows to estimate x from 
below, knowing only the spectral density g{E^) of one 
arbitrary eigenvalue E*, 

X = ^ g^iE) > g\E,) + ^ [1 - q{E.)\ = X ■ ( 3 ) 

E 

Here, we assume a completely flat density g{E) on its 
support aside from E*. x allows for a rather accurate 
extimation of x if is a single highly degenerate eigen¬ 
value compared to all other eigenvalues, i.e. if ^ 

g{E ^ E*), or if all eigenvalues are non-degenerate, i.e., 
p{E) = 1/iV for all E. These two limits correspond to 
vastly different networks. Eor instance, chain-like net¬ 
works with H belonging to the above mentioned class 
have eigenvalues with (mostly) the same degeneracy, i.e., 
p{E) = const, for all E, which yields X = X = 1/-^ [8], 
suggesting very efficient transport. In contrast to this, 
stars of the same size typically have a single highly de¬ 
generate eigenvalue yielding x = X = 1 ~ (4.^ — 6)/7V^ 
[8], thus, rendering transport inefficient. Obviously, in 
the infinite system size limit one has 


star-like topologies yield the same exponent k. This al¬ 
lows us to group networks, depending to their asymp¬ 
totic (global) quantum transport efficiency, into universal 
classes defined by k. 



FIG. 1. (Color online) Examples of three tree-like networks: 
(a) a scale-free tree (SFT), where the functionalities are drawn 
from the probability distribution P{fj = x) oc x~‘’ (here 
s = 2.5), (b) a Dendrimer (D) of generation 5 with func¬ 
tionality / = 3, and (c) a Viscek fractal (VF) of generation 3 
with functionality / = 4. Three different types of nodes are 
marked: leave nodes (green triangles), parents (brown circles) 
of leave nodes, and all other (open circles). 


^oo = ^ 

—oo Af-s-oo — 


0 for chains 

1 for stars 


( 4 ) 


Breakdown of Quantum Transport - For networks 
whose spectra dependent on a tunable paramter cr G M 
such that for large values of a the spectral density of the 
selected eigenvalue E* is of order 0{1/N) and for small 
values of a of order 0{1), one might observe a transition 
from efficient to inefficient transport. Let this transition 
occur at a given parameter value ac- At this value the 
transport breaks down which is reflected by a change of 
x(cr) from values of 0{1/N) for a > ac to values of 0{T) 
for a < Uc- This is reminiscent of a phase transition 
where, in our case, the quantity 1 — x (^) represents 
the order parameter. (We use the usual terminology of 
phase transitions in order to stress the similarities and 
to avoid to overload the paper with new terminology.) 
Consequently, we associate with this transition a critical 
exponent dehned by 


K = lim 
cr—>-crc 


log|l -X^(g)l 

log I CT - CTc I 


( 5 ) 


Since ^ Xoo(o'), one has k > n, where n is 

the exponent associated with the transition for Xoo(o'). 
As we will show below, tree-like networks which have 
a parametrized transition from chain-like topologies to 


Tree-like Networks - Figure 1 shows examples of three 
tree-like networks: (a) a scale-free tree (SFT), where the 
functionalities are drawn from the probability distribu¬ 
tion P{fj=x)(xx~^ (here s = 2.5), (b) a Dendrimer (D) 
of generation 5 with functionality / = 3, and (c) a Vis¬ 
cek fractal (VF) of generation 3 with functionality / = 4. 
For D, the functionality / refers to those nodes having 
more than one bond, and for VF, it refers to the building 
blocks of the fractal in each iterative step, see the squares 
as guide to the eye in Fig. 1(c). We do not distinguish 
between the respective functionalities, but we will always 
note by superscripts D and VF to which structure we re¬ 
fer to. D and VF allow for a direct computation of x and 
X based on Eq. (3), since their spectra are exactly known 
[10-12]. For SFT one has to resort to either numeric cal¬ 
culations of the spectra or to an analytic estimation of 
the spectral density of that eigenvalue with the highest 
degeneracy by counting those nodes (the leaves, green 
triangles in Fig. 1) of the network which only have a sin¬ 
gle bond. In doing so, we will also need to estimate the 
number of nodes (the parents, brown circles in Fig. 1) to 
which the leaves are connected, see below. 

As is easily shown, if a parent node j has two different 
leaves I and k, then the superposition state {\l) — |fc))/-\/2 
is a normalized eigenstate of the network. The total 
amount of these superposition states - all belonging to 
the same eigenvalue A* = H{fj = 1) = 1, regardless 
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their parents because we have assumed equal couplings in 
H - can be directly calculated. We note, that there could 
be also other eigenstates, not only involving leaves and 
parents, leading to the same eigenvalue [12, 13]. How¬ 
ever, since the number of these latter states is small we 
concentrate on the former states. Let be the num¬ 
ber of leaves of a single parent node j, then there are 
exactly — 1 independent eigenstates of such super¬ 
position type. Clearly, the more leaves a parent node has, 
the larger will be the number of independent eigenstates 
for a given eigenvalue, i.e., the larger will be the corre¬ 
sponding spectral density g{E^) > (A^l — Np)/N, where 
Nl is the number of all leaves and Np is the number of 
all parent nodes, which in turn can be related to N^: 
Assuming that a parent node j has functionality fj and 
that this node is also connected to Sj + 1 nodes which 
are not leaves, then the number of leaves of this node is 
{fj — 6j — 1). Let Ip be the index set of all parent nodes, 
then Np = J2jeir (/j “ “ 1)' ^ith Np = ^ 

we obtain 


Np 


Np 


( 6 ) 


where {f denote the av¬ 

erage over the set of all parents. Given the tree-like 
topology, we can estimate A^l based on N and the av¬ 
erage functionality of those nodes which are not leaves, 
{f)N\NL- For trees there are N — 1 bonds connect¬ 
ing the N nodes. In the two extreme cases of chains 
((/)Ar\jVL = 2) and stars ((/)jv\ArL = ^ - 1), one has 2 
and A^—1 leaves, respectively. For other trees the number 
of leaves 


Np = N- 


N -2 

if)N\Ni, ~ 1 


(7) 


which lies in the interval [2, A^ — Ij. We are now in the 
position of expressing q{E^) in terms of N and the aver¬ 
aged functionalities and {f)N\N^- Inserting into 

the RHS of Eq. (3), one obtains as a lower bound the 
final result up to order 1/A^: 


X > 



1 

(/ )n\Ni, ~ 7 



1 

(/-^)7Vp 
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1 - 


/ 2 \ 


/ (/-^)iVp-2 \ 


( 2 \ / ( / - - 2 
~ ~ '^)Arp ~ 7 


■ ( 8 ) 


In the limit when (/)jVp ^ i-e., when a 

parent node is only rarely coupled to more than one 
other parent node, one has — = (/)jvpj where 

(/)wp often be written as a function of (/)jv\ArL- 
If the functionality of non-leave nodes does not system¬ 
atically depend on the position in the network, one has 


{f)Np = {f)N\N^ = if)- We note, that there are ex¬ 
ceptions, see also the VF below. Equation (8) defines 
only a lower bound to x because we have neglected those 
eigenstates which are not simple superpositions of two 
states localized at leaves belonging to the same parent. 
In the infinite system size limit, these states are negligi¬ 
ble close to the transition point, such that the equality 
holds for Y . 

Considering the inverse of the average functionality as 
the network’s adjustable parameter, i.e., cr = !/(/), 
we can deduce a universal behavior at the breakdown 
of quantum transport. In the limit N ^ oo we obtain 

Xoo = (i-yTy^i) ’ (9) 


which results to first order inl/(/)inl — % =4/(/) 

and we get 


log[l-xJ 

l /</)^0 log[l/(/)] 


( 10 ) 


regardless of the original underlying network, be it de¬ 
terministic or random. Therefore, all tree-like networks 
will yield the same universal exponent k. We note, that 
the presence of loops could eventually lead to different 
exponents. 

Examples - In order to corroborate our general find¬ 
ings, we consider the three examples of tree-like net¬ 
works depicted in Fig. 1. All these networks allow 
for a parametrized transition from chain-like to star- 
like topologies, depending on the (average) functional¬ 
ity. While D and VF allow for a direct computation 
of X and x based on Eq. (3), we will employ Eq. (8) 
for SFT, which, depending on the average functionality, 
can have many leaves. For normalisation purpose it is 
inevitable for finite systems to impose a maximal func¬ 
tionality /max < V — 1 such that the average becomes 
(/) = Inserting this average in 

Eq. (8) yields We note that (/) depends on the 

scaling parameter s and that !/(/) —>■ 0 when s \ 2. 
In limit V —>■ oo the averages are related to the Rie- 
mann zeta function /(s) = allowing to write the 

leading terms for values of s > 2 as 


X®^^ = l-4^ 


C(s) - 1 


C(s- 1) -C(s) ' 

Therefore, the critical exponent follows as 
log(l - x®’"'^) 

K = lim -;-;-—- = 1 , 

s\,2 log(s - 2) 


( 11 ) 


( 12 ) 


which again confirms our general statement about the 
universal behavior of tree-like networks. Figure 2 shows 
the dependence of 1 — x®^"*" iF® scale-free parameter 
s and, as inset, on (/) for different sizes N and also 
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for N ^ oo. For finite SFT, we have compared our 
analytic estimation given by Eq. (8) (lines) with Monte- 
Carlo simulations (symbols) for SFT grown iteratively by 
the algorithm given in [14] with the connectivity matrix 
defining H. In the numerical computations we have con¬ 
sidered ensemble averages of with ensembles sizes of 
R= 10®/fV. All curves show the expected scaling close to 
the transition point where 1/ ( /) —t 0 and s \ 2, respec¬ 
tively. One notes from the inset of Fig. 2 that with in¬ 
creasing N the finite-size effects become less pronounced, 
leading eventually to a sharp transition for s 2. 



FIG. 2. (Color online) Numerical (symbols) and analytical 
(lines) results for 1 — For finite N = 10^,10^, and 

10^, the ensemble averages for 1 — with the connectivity 
matrix defining H and with R—IO^/N realizations are shown 
(in linear scale) as functions of the scale free parameter s. The 
analytical estimates are obtained from Eq. (8) with <5 = 0 for 
finite SFT and, in the limit N —^ oo, from Eq. (11). The inset 
shows the same curves (in log-log scale) as functions of 1/ {/). 
In both plots, one notices the clear sign of the breakdown 
of the quantum transport, with the expected linear behavior 
close to the transition at s \ 2 or, equivalently, 1/ (/) —>-0. 


Eq. (13), we obtain for 1/ ( / ) —0 that 1 —X° ~ 1/ if) 
and that 1 — x^^ ~ 1/ if) ■ Thus, also here we find the 
breakdown leading to the (exact) exponent k = k= 1. 

We finally want to stress the differences between D and 
VF on the one hand and SFT on the other hand: For 
all structures, the transition happens when 1// —>■ 0 or 
1/ (/) —0. However, for D and VF, having fixed deter¬ 
ministic functionalities /, there is no parameter allowing 
to study the behavior of Xoo beyond the critical point. 
Moreover, the limit / —> oo seems rather artificial in the 
TV —>■ oo limit since no real system can ever reach this. 
The situation is different for SFT: Even though (/) di¬ 
verges for scaling parameters s < 2, it is possible to study 
the behavior of in fhi® parameter region. We fur¬ 

ther note that for finite SFT, one observes the maximal 
values of < 1 only in the limit when s —>■ 0. 

Conclusions - We have shown, that the breakdown of 
quantum transport on complex tree-like networks shows a 
universal behavior characterized by a global transport ef¬ 
ficiency measure based on the time-averaged return prob¬ 
ability. Parametrizing the corresponding Hamiltonian by 
the (average) functionality of the nodes of the network 
allows to derive bounds for this measure, leading in the 
infinite system size limit to a characteristic universal ex¬ 
ponent for all tree-like networks. We anticipate that a 
similar treetment might also be feasible for other net¬ 
works, also including loops. 
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For D and VF, the exact knowledge of their spectral 
densities [10-12] allows us to calculate x as well as its 
lower bound x as functions of the respective functional¬ 
ities. According to Eq. (3), we get in the infinite size 
limit 


X^= [I-j 


and x^^ = 1-6- 


/-I 


/(/ + 2)-2 
as well as the corresponding lower bounds 


(13) 


X^= 1- 


1 


/-I 


J VF I 1 4/ 5 

and x^ = (1 - 


(14) 


In both cases we find the breakdown of transport in 
the limit 1// —0. In order to be comparable to the 
SFT, we express both, xS> and x^^j as functions of ( /). 
For D, ( / “ /’ which is not the case for VF, where 

if)N\N^ = (/ + 4)/3 while (/)^^ = /. Inserting into 
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